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We present a detailed quantitative comparison between theoretical calculations of the local density 
of states and recent experimental measurements of scanning tunneling spectra around Ni impurities 
in Bi2Sr2CaCu208+a . A double-peak structure on the hole side of the spectrum at the Ni site 
is identified as the spin-split van Hove singularity in the band structure. The Ni atom induces 
local changes in hopping matrix elements comparable in size to the induced on-site spin-dependent 
potential. We find evidence from the measurements of order parameter suppression in the vicinity 
of the Ni impurity. These extracted impurity parameters can be of use in quantitative calculations 
of macroscopic response properties, such as the AC conductivity. 



PACS numbers: 74.25.Jb, 74.80.-g 

The response of a superconductor to impurities can 
determine several macroscopic electromagnetic proper- 
ties, such as thepjnagnetic penetration depth and high 
frequency losses.EI For high-temperature superconduc- 
tors it is advantageous to study the impurity proh-> 
lem in real space rather than in momentum space,cl 
not only because the impurities destroy the transla- 
tional invariance, but also because these doped systems 
are intrinsically inhomogeneous.tj High-quality electronic 
spectra, using the technique of scanning tunneling mi- 
croscopy (STM), around individual impurity atoms in su- 
perconducting I^S^CaC^Os+a have recently become 
available. Quantitative comparison between calcula- 
tions of the local density of states (LDOS) and these STM 
measurements can lead to effective impurity parameters 
that are useful for calculating macroscopic properties. 

A number of theoretical studies of the STM spec- 
tra near impurities in d-wave supercondujE.toga-have ,al-. 
ready been carried out by various authors.l3i3'l0&BE3O-E£l 
These previous studies focused on qualitative features, 
such as the existence of quasiparticle resonances, the 
spatial structure of the resonances, and the differences 
among systems with various pairing symmetries. Often 
the model for an impurity only consisted of effective po- 
tentials at the impurity site, or relied on an unrealistic 
band structure for the host superconductor. There are 
still significant quantitative disagreements between the 
theoretical and experimental LDOSJH3'EJE3 Of the types 
of impurities and defects for which high-quality STM 
spectra are available, the Ni impurity provides a good 
opportunity to resolve some of these issues. Because the 
perturbation to the-, local environment caused by a Ni 
impurity is weak,Elt2l we find that, with some extensions 
to the potential model the STM spectra can be well re- 
produced by our calculations. In addition to the short- 
ranged potentials at the impurity site, we also include 
the change of the hopping matrix elements from the im- 
purity site to its nearest neighbors, and the suppression 
of the order parameter. 

We also find an interesting interplay between the un- 
derlying band structure and the superconducting coher- 



ence peak near a Ni impurity. In the spectrum of the Ni 
site there is a large peak on the hole side (below the chem- 
ical potential), right outside the gap edge. We find that 
this peak, previously thought as the remnant of the co- 
herence peak, is largely contributed by a density-of-states 
peak in the underlying band structure. One special and 
important feature of hole doped cuprate superconductors 
is that a van Hove singularity is close toJi^F^piiiJ^el. 
It has been suggested by various authorsBBBBBlia'El 
that the existence of the van Hove singularity could affect 
the mechanism of superconductivity, and also account 
for the unusual normal state properties. The existence 
of such a density-of-states peak has received direct sup- 
port from experiments using .angle-resolved photoemis- 
sion spectroscopy (ARPES).clrE3 Despite this density-of- 
states peak (aside from the coherence peaks) in the super- 
conducting state has never been directly reported in tun- 
neling experiments. Traces of this feature were sugg ested. 
to explain the asymmetry of the coherence peaksJHj'Eil 
and the sharpness of the coherence peaks. However, no 
specific qualitative influence of a van Hove peak on the 
midgap impurity_t|esonances of a particular impurity has 
been considered .l3 

We find that the Ni impurity offers a novel probe of the 
van Hove peak. The van Hove singularity is likely very 
close to the Fermi level, therefore, in the superconducting 
state the van Hove peak cannot be resolved from the 
coherence peak. The Ni impurity will locally suppress the 
coherence peak, but only weakly perturb the underlying 
band structure. As a result, the van Hove peak shows up 
unmasked in the spectrum at the impurity site. 

Our calculation is based on the following tight-binding 
Hamiltonian, 

H = - + St v) c L c i<r 

+ (V + Vs)4 t cdt + (Vo - V s )cl lCol , (1) 
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FIG. 1: A schematic diagram showing our model of a Ni im- 
purity on the Cu-O plane. Parameters extracted from fitting 
to the data are listed on the right. 



where i and j label sites (the impurity resides at site 
0), and a labels spin. The hopping matrix elements, 
Uj, are taken—from a one-band parameterization of the 
ARPES data.EEl The numerical values in units of meV are 
148.8, -40.9, 13, 14, and -12.8 for the hopping to the 
first five nearest neighbors. The measured value for the 
gap maximum Ao at a site away from the impurity is 28 
meV.B The order parameters, A^ , are assumed to be real, 
and to have d-wave symmetry. They are only non-zero 
on the bonds between two nearest-neighbor sites, and 
A^i+x = — Aj^+g = Ao/4. The momentum-dependent 
order parameter resulting from a homogeneous system 
with these Ay is Ak = (Ao/2)(cosfc :r — cosk y ), where 
the lattice spacing between two Cu atoms is unity. 

Because the experiments were carried out at low tem- 
peratures, the spin degree-of-freedom of the magnetic Ni 
impurity is assumed to be completely frozen. At the Ni 
site, we use a potential that consists of a magnetic part, 
Vs, and a non-magnetic part, Vq. It is found that the 
Ni impurity induces a change in the hopping, dUj , to the 
nearest neighbor sites, and a suppression of the order pa- 
rameter magnitudes, <5Ay , on the four bonds connected 
to the Ni site, and on the other twelve bonds connected 
to the nearest neighbor sites. A schematic diagram of our 
model and the fitting parameters used for the Ni impurity 
are shown in FIG. [l]. 

Calculations presented here are based on a Koster- 
Slater techniquc.EjEZlnfj The homogeneous Green's func- 
tion, g(rj, Tj, uj), in the Nambu formalism is first con- 
structed with a 1 meV quasiparticle linewidth. In 
the presence of an impurity modeled by V, the 
Gor'kov equation, G = g + gVG = (I — gV) _1 g, 
is solved by numerically inverting the matrices. The 
LDOS is the imaginary part of the Green's function, 
(-l/7r)Im[Gn(r,,r,,w) - G* 22 {r l: r,, -uj)]. 

In the following, we will discuss how various parame- 
ters influence our calculations. Although six parameters 
including the chemical potential are introduced in the 
model, there are particular features in the LDOS spectra 
associated with each of them. This allows us to determine 
each parameter largely independently from the others. 
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FIG. 2: LDOS spectra at various sites near a Ni impurity. 
The left panels show the spectra in the normal state. The 
right panels show the spectra in the superconducting state. 
Solid lines show the calculated results. Dashed and dotted 
lines show the two spin components, (a) and (b) are the 
spectra at a site about 30A away from the Ni impurity, (c) 
and (d) are the spectra right at the impurity site, (e) and (f) 
are the spectra on the nearest neighbor site, (g) and (h) are 
the spectra on the second nearest neighbor site. Open circles 
(o) show the STM data.S The data was rescaled by a constant 
factor identical for all the spectra. 



Let us first consider the.-Ghemical potential. From the 
ARPES parameterizationpj the chemical potential for 
an optimally doped sample is —130.5 meV. However, we 
know from STM measurements that the doping concen- 
tration is locally varying, andjs correlated with the su- 
perconducting gap magnitude.0 A region with a gap mag- 
nitude of 28 meV is overdoped. The fit to STM spectra 
is optimal if the shift in chemical potential relative to the 
optimally doped region is about —25 meV. This amount 
of shift changes the location of the van Hove singularity 
from 34 meV to 9 mei/J below the Fermi level, shown in 
FIG ^(a) for Ay = 0.c3 In the superconducting state the 
van Hove singularity cannot be resolved from the coher- 
ence peak at the hole side, shown in FIG ||(b). The chem- 
ical potential shift also sharpens the coherence peaks, and 
contributes to the asymmetry between the two peaks. Al- 
though not direct evidence for the van Hove singularity 
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FIG. 3: Evolution of the spectra as a function of distance 
from the Ni atom. Left panel shows the calculations, and the 
right panel shows the data from Fig. 4 in Ref. 5. The spin 
splitting of the resonance peaks inside the gap is constant as 
one moves away from the impurity site. The spin splitting of 
the van Hove singularity decreases as the distance increases. 



in STM, the agreement between the calculations and the 
measurements is certainly improved. 

We now describe the modeling of the Ni impurity fo- 
cusing first on the single-site potential model with only 
the on-site potentials, Vq and Vs- For a weak potential, 
the main effect of these potential is to set the energies of 
the resonances. Since the resonance peaks for Ni are at 
9.2 and 18.6 meV, and the perturbation from Ni is weak, 
both Vb + Vg and Vq — Vs should be negative. Our Vq dops. 
not differ greatly from that found in a simple modcl.Ej 
The two on-site parameters, Vq and Vs, are both required 
to fit the energies of the two resonance peaks.Eil Note that 
we have no control over the peak height or width in the 
spectra using this simple model. We find, due to ex- 
tended states within the gap, that the resonance peaks 
will be broader than the coherence peaks. In the observa- 
tion, however, the width of the resonance peaks appears 
comparable, or even sharper, than the coherence peaks. 

In order to reduce the coupling to the extended states 
and increase the localization of the resonance states, 
we then assume that the hopping to the nearest neigh- 
bor sites, 6tox, is reduced. With only a 20% change 
in the nearest-neighbor hopping, we can reproduce the 
linewidth of the resonance peaks. Vq and Vs must be ad- 



justed slightly when Stoi is introduced in order to keep 
the resonances at the fixed energies. Our calculations 
show that the magnitudes of the on-site potentials change 
to even weaker values, and characterize an even weaker 
impurity potential. Using these three parameters for cal- 
culating the effect of a Ni impurity in the normal state 
(shown in FIG ||(c)), one can see that the band struc- 
ture is not highly distorted. The van Hove peak remains 
reasonably sharp, and splits into two spin components. 
Comparing ^(c) to the superconducting state, shown in 
FIG ||(d), we can identify the two peaks on the hole side, 
right outside the gap edge, as the spin-split van Hove 
peaks. A similar connection between the normal state 
electronic structure and an impurity's influence on the 
superconducting state has been identified ip-caJculations 
for impurities in s-wave superconductors .E3c3 Another 
distinction between the midgap resonances and the van 
Hove peaks is that the spin splitting of the resonance 
peaks inside the gap remains independent of distance 
from the impurity site. The spin splitting of the van Hove 
singularity, however, decreases as the distance increases. 

Once the large peak at the hole side was identified as 
the van Hove peak, a noticeable disagreement between 
the calculation and the data became apparent. The cal- 
culated relative energy difference between the resonance 
peaks and the van Hove peaks is somewhat larger than 
the observed value if the gap maximum is constant every- 
where. In the d-wave superconducting state, the position 
of the van Hove peak is dominated by the gap edge. One 
cannot increase the agreement by tuning Vq because that 
shifts all peaks almost an equal amount of energy. We 
could obtain good agreement if the gap maximum were 
about 22 meV instead of 28 meV. In principle using self- 
consistent calculations can account for the effect, how- 
ever, the self-consistency conditions involve the pairing 
mechanism, which is controversial. Taking a phenomeno- 
logical approach, we can identify the suppression of the 
order parameter around the impurity from the data. This 
then becomes a constraint on pairing theory. We find 
that the parameter is locally suppressed more than 50% 
(See FIG. 0). 

So far, we were only fitting the spectrum right at the Ni 
site. Nevertheless, good agreement in spatial structure is 
also obtained. Detailed comparisons for the spectra on 
neighboring sites are shown in FIG. ||(e)-(h), the spatial 
decay of the resonance peaks is shown in FIG. 0, and 
the spatial maps at given energies are shown in FIG. |J. 
These are direct comparisons with the data in FIGs. 2-4 
in Ref. 5. 

Two key additional features are found from the spa- 
tial structure. The first one is that the.iiple component 
of the resonance has d-wave symmetry.BEJ The second is 
that the perturbation, due to the van Hove singularity, 
is strongest along the diagonals of the lattice. At given 
energy, the secondary peaks are always located at the 
next-nearest rather than the nearest neighbor sites rela- 
tive to the main peaks. This is because, around the Fermi 
level, the saddle points that give rise to the van Hove sin- 



4 




FIG. 4: (color) Spatial structure of the LDOS around a 
Ni impurity in the normal state at (a) —15 meV, (b) —50 
meV (near van Hove peak positions in FIG. [](c)), and in the 
superconducting state at (c) +9 meV, and (d) —9 meV. The 
Wannier function is assumed to be a Gaussian with a width 
of half the lattice spacing. The LDOS is shown in logarithmic 
scales and in units of eV _1 . The lattice orientation in (c) and 
(d) needs to be rotated by 45° when compared to Fig. 2 in 
Ref. 5. 



gularity in the density of states are near the (±tt, 0) and 



(0, ±7r) points in momentum space. The main contribu- 
tion to the real space Green's function largely comes from 
these high density-of-states regions, 



(2tt)2 



= jk-(r; 



g(k,w) 



(2) 



For Yi — Tj = (±1,±1) the four regions add up con- 
structively, while for (±1, 0) or (0, ±1), the contributions 
from (±7T, 0) cancel with those from (0, ±7r). As a result, 
|g((l, 1), (0, 0), w)| » |g((l, 0), (0, 0),u)\ for energies near 
the chemical potential. This is also consistent with the 
larger spin splitting of the van Hove peak at the second 
nearest neighbors compared to the nearest neighbors. 

In this study we have assumed that the STM spectra 
are directly related to the electronic structure of the Cu- 
O plane regardless of the presence of the intermediate 
Bi-0 plane. As was argued in Ref. 5, the energies and 
linewidths of the resonances should be the same. We find 
the principal effect of a d-wave filter^ on our parameter- 
ization is to (roughly) change the sign of Vq. Therefore, 
the evidence for the reduction of hopping, and suppres- 
sion of the order parameter are still clear. 

In summary, the presence of a Ni atom on the Cu- 
O plane produces a weak perturbation to the system. 
We show direct evidence for a density-of-states peak in 
the STM measurements and local suppression of the or- 
der parameter. The suppression of the order parameter 
can constrain the self-consistency condition and thus the 
pairing mechanism. 
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ing data shown in FIGs. || || from Ref. 5. This work is 
supported by ONR Grant No. N00014-99-0313. 
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